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Glass-melting processes and the effect of batch quality on these processes are considered. An analysis is per¬ 
formed on the basis of a heat-engineering study of large-capacity glass-melting furnaces and mathematical 
modeling. Admissible single-impulse and multi-impulse deviations from a prescribed batch formula are esti¬ 
mated. 


Highly efficient float-glass production lines and new ma¬ 
chines for molding hollow articles that have replaced old 
glass-molding machines require substantial intensification of 
the performance of glass-melting furnaces. The need for in¬ 
crease furnace efficiency is also motivated by economic fac¬ 
tors (growing prices of refractories, fuel, and power) and 
market factors (growing competition of manufacturers, 
stricter requirements imposed on product quality). The most 
essential factor in the production of sheet glass is the fact that 
it has become a basic product for making multiple window 
panes, laminated building and car glasses, and coated glasses 
(mirrors, energy-saving, decorative glass), i.e., articles that 
have to meet special requirements. 

The use of float glass as a finished product has steadily 
decreased. In 2002 Russia produced 5.5 million m 2 of glass 
with coatings, 2.2 — mirrors, 6.3 — laminated glass, 17.0 — 
multiple widow glasses, and 8.3 million m 2 hardened glass, 
which totals 65% of basic float glass produced. 

The efficiency of furnaces for container and sheet glass 
has reached 2000 kg/m 2 of the working zone of the furnace 
and even more. The furnace campaign lasts 7 years and longer. 

In these conditions, requirements imposed on raw mate¬ 
rials and batches have grown perceptibly. The low quality of 
batches and initial materials is responsible for the fact that 
production of high-quality sheet and hollow glass in Russia 
is limited and substantial quantities of these products and 
materials are imported, which presents challenging competi¬ 
tion to domestic manufacturers. 

Let us analyze glass-melting processes and the effect of 
batch quality on glass melting based on some results pub¬ 
lished earlier [1], 

Figure 1 presents a scheme of batch melting in a tank 
furnace. 

Batch is charged into the furnace on top of glass melt at a 
velocity v (m/h). As the batch layer keeps moving along the 
furnace, its thickness changes due to its melting from above 
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and from beneath. The melt emerging on the top is running 
down from the tops of batch heaps, whereas the melt emerg¬ 
ing beneath directly enters the glass melt. Within 1-2 min 
after batch loading, a viscous melt film of thickness 10 mm 
is formed on the surface of the batch, which impedes its com¬ 
plete melting. It can be assumed that the process of heating 
and melting of the batch proceeds in a quasistationary mode, 
when the melt formation rate is equal to the rate of its run¬ 
ning down the batch heaps. The same assumption can be ex¬ 
tended to the process of batch heating from beneath. 

Let x and y be immobile coordinates. Then the following 
dependences are true for silica dissolving in a silicate melt 
(silicate and glass formation): 
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Fig. 1. Batch melting scheme. 
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t m m 

the time of charging start; J Wy 0x is the quantity of 

t 0 J=' i=1 

Si0 2 fixed as a result of silicate-formation reactions; m is the 
number of components; K ; is the constant of the zth reaction 
of silicate formation. 

Equation (1) can be complemented by the heat exchange 
equation: 

pc' — = div (A grad T— qDp grad w .) - pv grad H, 


where c ’ is the heat capacity of the mixture (silicates + 
quartz + melt); c' = (p Sl02 w Si02 + <p m w m (<P si02 and <p m are 

the contents of Si0 7 and melt in the silicates + quartz mix¬ 
ture); A is the effective thermal conductivity of the melt; q is 
the specific isothermic transfer heat; w m is the weight content 
of the melt in the mixture with dissolved Si0 2 particles; H is 
the enthalpy of melt mixture, H = H w si0 ^ + H m w m . 

Equation (1) was analyzed using data obtained in a heat¬ 
engineering study of a float-line furnace (Saratov Technical 
Glass Works) and a vertical drawing line (Mekhsteklo Lisi- 
chanskii Works). 

The rate of change of Si0 7 content in the longitudinal 
section of the furnace is approximated by the dependence 

w Si0i = 13.75 + 4.5x + 0.14y- 1.2/- 7.8x 2 + 1.6/ 2 , 

v a 8 - 15 m/h, D a 10“ 7 cm 2 /sec. 

Then the first term on the right-hand side of Eq. (1) has 

<3wo in „ 

the order v- -~40-80; v-150-250, the 
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second term has 10 - 4 , and the fourth term 10 - 30. 

The estimate given above makes it possible to simplify 
the problem to a unidimensional case and to obtain the fol¬ 
lowing transcendent equation with respect to 



where x = S is the coordinate of the glass melt - friable batch 
boundary: 
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where j is the specific heat of batch melting. 

The following data were assumed for the calculation of 
the melting zone boundary: 

- the batch-charging velocity for the float-line furnace is 
16 m/h and for the vertical-drawing line 10 m/h. 

- for the batch composition corresponding to the glass 
formula: A b = 0.1914 W/(m • K), A g = 132 W/(m • K), and 
j = 8.36 J/kg. 

- the batch layer thickness (2/3 )L = 0.2 m, T m = 1123 K. 


Then based on Eq. (2) we obtain that the duration of 
batch heating up to the melting point is approximately 
0.67 h. Considering the batch charging velocity in the Lisi- 
chanskii Works furnace, this corresponds to the melting zone 
boundary of 6.7 m (the 2nd burner zone), for the Saratovskii 
Works furnace it is 11.4 m (between the 3rd and 4th burners). 

The initial parameters and, consequently, calculation re¬ 
sults depend on the melting temperature and stability of the 
batch composition. Under a constant melting temperature, a 
2% deviation in Si0 2 content changes the calculated melting 
zone boundary by ± 1.3 m. 

It is known that a steady position of the melting zone 
boundary is one of the most important factors for maintain¬ 
ing a stable melting regime and ensuring glass melt quality. 

The stability of the chemical composition of the batch, 
other terms being equal, depends on the accuracy of compo¬ 
nent proportioning. The proportioning accuracy, in turn, de¬ 
pends on the proportioner system design and the quality of 
the electronic control system. Furthermore, a preset batch 
composition may differ from an actual composition due to 
the fact that the materials (batch components), as a rule, in¬ 
troduce several glass components (oxides) simultaneously. 

It is possible to define more strictly the following tech¬ 
nological and physical limitations in estimating the propor¬ 
tioning accuracy [2-4]: 

- rating of glass components; 

- some glass components are introduced by means of 
several batch material components; 

- in the system of equations describing the dependence 
of the batch composition on a preset glass composition, the 
quantity of unknown values may differ from the number of 
equations (the case of a nonquadratic matrix); 

- the solution of the system of glass composition - batch 
composition equations should be non-negative; 

- it is necessary to take into account mechanical entrain¬ 
ment of batch components and volatilization of glass compo¬ 
nents from melted glass; 

- compositions of materials, batch, and glass are influ¬ 
enced by errors depending on measurement methods and in¬ 
strumental errors. 

The material composition - batch composition depen¬ 
dence is described by the following relationship: 

Y = AX, 

where Y= {y .} is the //-dimensional vector of the batch 
composition; X= {xj } m is the /n-dimensional vector of mate¬ 
rial portions; A = {a ?/ } is the matrix of dimensionality of the 
material composition n x m. 

Passing on to a differential form and rejecting 2nd order 
terms, we will obtain 

°NZ[(Vi/) 2 + (V/) 2 ]’ (3) 
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where a i is the mean quadratic deviation Ay.., AjA- = 
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Fig. 3. Dynamics of admissible changes in the batch composition. 


Az, % 



It can be seen that only the DSI proportioner ensures pro¬ 
portioning accuracy meeting the requirements of the Techni¬ 
cal regulations for sheet and container glass factories. 

A single impulse change in a batch composition was con¬ 
sidered above. For a multi-impulse deviation in Si0 2 content 
in a batch, another approach to estimating an admissible de¬ 
viation level is possible [5]. 

Figure 2 shows fluctuations in SiO, content in a batch 
and respective fluctuations in Si0 2 content in glass produced 
at the Tokmakskii Works: 


Fig. 2. Variations in Si0 2 content in the batch Aw within time pe¬ 
riod At and respective variations in SiO, content in glass Az: dashed 
curves) variations in SiO, content in glass for a single impulse 
change in Si0 2 content in the batch; solid curves) variations in SiO, 
content in glass for multi-impulse change in SiO, content in the 
batch. 
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tent of the /th material component; o^ = ^Ax i2 a j (a j is the 

accuracy of detennining the content of the main oxides in the 
y'th material component); Ax - is the absolute error of mea¬ 


suring the content of the ith oxide in the /th batch compo¬ 
nent; A A • is the absolute error of measuring the weight of 


theyth material component; x ; is the measured content of the 

1 - 

ith oxide in the /th material component; ct = - A A • max . 


Calculation based on relationship (3) was carried out for 
two types of proportioners used in the Russian glass indus¬ 
try: a DVST proportioner (with a lever system and weight 
fixation based on the indicator pointer) designed by the Veda 
Company and a DSI strain-measuring gauge (with a MBUD 
microprocessor controlling block) designed by the Institute 
of Glass JSC and the Stromizmeritel’ Company (Table 1). 


TABLE 1 


Weighing error, % of the maximum weighing limit 


tioner 

sand 

dolomite 

chalk 

sulfate 

soda 

feldspar 

DVST 

0.26 

0.66 

0.83 

0.67 

0.32 

1.06 

DSI 

0.12 

0.16 

0.31 

0.28 

0.14 

0.32 


A adm.b 


A adm.g 


£a/ a 2 , (iAl) 


where A admb and A admg are admissible SiO, fluctuations in 
the batch and in glass; A f Al is the change in the SiO, content 
in glass (a result of modeling the homogenization process); 
i is the number of impulses. 

The calculation results are indicated in Fig. 3. 

Thus, a single deviation (one batch portion) up to 1% can 
be permitted due to the homogenizing capacity of the fur¬ 
nace. As a rule, the deviation should not exceed 0.15%. Such 
weighing accuracy can be provided only by the DSI propor¬ 
tioner or imported dosing machines with a class of accuracy 
not lower than that of the DSI. 
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